Non-local order parameters for the ID Hubbard model 
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We characterize the Mott insulator and Luther-Emery phases of the ID Hubbard model through 
correlators that measure the parity of spin and charge strings along the chain. These non-local 
quantities order in the corresponding gapped phases and vanish at the critical point Uc = 0. The 
Mott insulator consists of bound doublon-holon pairs, which in the Luther-Emery phase turn into 
electron pairs with opposite spins, both unbinding at Uc- The behavior of the parity correlators can 
be captured by an effective free spinless fermion model. 
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The Hubbard model and its extensions have been 
widely used to investigate the behavior of strongly corre- 
lated electrons in several condensed matter systems rang- 
ing from Mott insulators (MI) to high- Tc superconduct- 
ing materials. Recently, the progress in ultracold gas ex- 
periments that use fermionic atoms trapped into optical 
lattices has opened the way to the direct simulation of 
the Hubbard model and the observation of the predicted 
MI phase [T]. Since the Mott transition is of Berezin- 
skii Kosterlitz Thouless (BKT) type, the MI phase does 
not admit a local order parameter; instead the transition 
point corresponds to the vanishing of some topological 
order, possibly described by appropriate nonlocal quan- 
tities [H|3j. A progress in this direction has been achieved 
in the related field of the bosonic Hubbard models, where 
MI and Haldane insulator phases have been character- 
ized by means of non-local string parameters, inspired 
by the correspondence of the bosonic system with spin- 
1 Hamiltonians at low energy near integer filling jJl [S] . 
One of these parameters is related to the parity correlator 

Op{r) = (e^"^J=' '^"■'), with S^^i = \{ni-v) measuring 
the parity of the deviation of the occupation number n,j 
with respect to the filling in a string starting from the 
site i, ending to the site i + r. The non-vanishing value 
of the parity parameter Op — lim^-^oo Op{r) in the in- 
sulating phase has been observed with in situ imaging in 
experiments on ultracold bosonic *^Rb atoms [6 . 

In this Letter we address the study of nonlocal string- 
type correlators to inspect the gapped phases of the 
fermionic Hubbard model. The expected role of antiferro- 
magnetic (AF) correlations has so far driven the attention 
mainly to the study of Haldane type string correlators; 
these were found to vanish algebraically, together with 
Op{r) in the Luttinger liquid regime [7]. On the other 
hand, in the large Coulomb repulsion limit the Hubbard 
Hamiltonian at half-filling is known to reduce to the AF 
Heisenberg Hamiltonian, for which the parity string cor- 
relator reduces trivially to the identity, the wavefunc- 
tion being frozen to the sector with only one electron per 
site. Since in the MI phase the number of doubly oc- 
cupied sites (doublons) and empty sites (holons) is non- 



vanishing at any finite value of the interaction (as also 
observed experimentally [T]), it is reasonable to expect 
that an appropriate parity parameter could characterize 
the crossover from the Heisenberg to the Luttinger liquid 
limit, marking the existence of the MI phase. 

The local 4-dimensional vector space on which an elec- 
tron Hamiltonian acts is typically generated by applying 
to the vacuum operators forming a su(4) algebra, with 
three Cartan generators. Consequently, we can introduce 
two independent parity correlators Op \ defined as: 
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with index v — c, s, namely the "charge" and "spin" gen- 
eralizations of the parity correlator Op{r). Here S^'^^ 
are the spin and pseudospin operators defined respec- 
tively as S^!"} = jl'^i.t ~ '^i.i) ^i^J = - l)i with 
n-ia — cl^Cia, cr ~tii: c^cr Creating a fermion at site i with 
spin CT. By means of bosonization and DMRG analysis, 
we will show that each Op orders in the corresponding 
gapped phase: MI for v — c, with open charge gap, and 
Luther Emery (LE) for i> — s, with open spin gap. The 
Op'' vanish with the gap at the BKT transition point 
where the correlation length becomes infinite. 

The Hubbard model is described by the Hamiltonian 



H = - 



(2) 



where the overlap integral U gives the on-site contribu- 
tion of Coulomb repulsion, and energy is expressed in 
units of the tunneling amplitude. 

The bosonized form of the half-filled Hubbard Hamilto- 
nian at low-energy is known to give rise to two continuum 
models describing separately the spin and charge sectors 
[8J. The latter is described by the Hamiltonian 



Hr = dx 
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with 

«c - VF (l + —) ' K,= (l + —) ' . (4) 

Here $c is the compactified boson describing the charge 
excitations with velocity Vc, and lie — dxQc/'^ is its con- 
jugate momentum (a is a cutoff). At the BKT transition 
point [/ = 0, we have K^, = f. The bosonic field in the 
spin sector $s is governed by equations which can be ob- 
tained from (|3| and Q by replacing U ^ —U and c ^ s. 
The spin-charge transformation c^^ — > (— l)-'cj^, that im- 
plies S'f^i — S'^'*], in the present bosonization analysis 
corresponds simply to the change <i>c $s. In fact, we 
have used the continuum prescriptions used in Ref. [8 
where ^-(x) = and J^(x) = 

v27r yz-K 

For [/ > 0, we get Ks > 1: the cosine term in Hg 
is (marginally) irrelevant and the spin excitations are 
gapless and governed by an ordinary Gaussian model. 
Meanwhile, Kc < 1 and a charge gap is generated by 
the relevant cosine term in He- As a consequence, the 
field <I>c is pinned in one of the classical minima of the 
cosine term, i.e. $c = G Z, while $s does not 

order. For U < 0, just the same occurs with inverted 
roles <i>c ^ $s- In the continuum limit one can realize 
that the parity operators become [U |9] 

oPir) « (cos[y2$,(r)]cos[y2$,(0)]). 

Hence in the MI phase ai U > 0, Op turns out to be 
non vanishing. In the U < case instead the LE phase 

(s) 

is characterized by nonzero Op . The two Haldane type 
string correlators O^ {r) = (^"..e'^^^S 

give instead OP{r) w (sin[v/2$^(r)] sin[\/2$i.(0)]) 
where the same argument suggests that these are both 
asymptotically vanishing in the two gapped phases. From 
the above derivation, we can conjecture that a necessary 
and sufficient condition for having an asymptotically non 
vanishing charge (spin) parity correlator in the Hubbard 
model is the opening of a gap in the charge (spin) sec- 
tor, so that Op ■* do configure as order parameters for the 
gapped phases of the Hubbard model. 

Below we support our previous argument providing a 
quantitative estimation of the parity string parameter in 
the MI phase. This is achieved by means of numerical 
analysis using the density matrix renormalization group 
(DMRG) algorithm on finite size chains with periodic 
boundary conditions (PBC's). The analysis requires very 
precise and reliable data; in fact, the computing effort is 
significant due to both the slowdown caused by PBC's 
and the high sensitivity of the correlations contained in 
Op''(r) with respect to numerical errors. Hence we have 
chosen to consider chain sizes from i = 10 to i = 50 
and 1024 DMRG states. The curves of 0''p\r) plotted 
in Figjl] for L — 50 evidence clearly a fast convergence 
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Figure 1: Parity correlator Op{r) for a periodic chain with 
1/ = 50 as a function of the string length r. The sequences 
of data refer to U = 0.1,1.0, 2.0, 3.0, 5.0, 10.0 (in ascending 
order). 

to the asymptotic values for high interactions as well as 
a progressive increase of the parity order with U. The 
presence of two sequences for even and odd r that tend 
toward the same asymptotic limit also signals that the 
spin parity correlator Op\r) — (— l)''Op''(r) has a uni- 
form part [O'^p^ (2r + + Op^ (2r)]/2 that goes smoothly 
to zero for U > 0. The opposite mechanism holds for 
negative values of the interaction. 

Exactly at [/ = both parity orders are absent and 
Op\r) — Op\r) as required by the spin-charge sym- 
metry. Here, an analytic calculation of 0^p\r) can be 
performed independently for both spin species by us- 
ing the Wick theorem and evaluating Toeplitz determi- 
nants. An estimation of the asymptotic behavior gives 
OP{r)r^r~^ a.tU = Oimi- 

(c) 

We have explicitly evaluated the order parameter Op 
in the MI phase and plotted it in Fig[2] for several val- 
ues of U. The asymptotic values have been extrapolated 
from the finite-size scaling of the quantity Op\L/2) in 
a periodic chain of length L. For the fits, we have made 
use of functions Op{r) = Op + Ar ^e^^/^ obtaining a 
good convergence. Interestingly, as evidenced in the in- 
set of Figj2j for small U we get 7=1 and A > 0, and 
for strong interactions we obtain 7 = 1/2 and A < 
while for intermediate values the best fit seems to be a 
combination of the two functions. 

The non vanishing of Op implies the existence of 
bound doublon-holon pairs; their correlation length in- 
creases by decreasing U becoming infinite at the tran- 
sition, when pairs finally unbind. The quasi long-range 
AF order of the MI phase suggests that such pairs are 
diluted in an AF background of single electrons. The 
spin-charge transformation that maps positive U Hamil- 
tonian at half-filling into negative U case at zero mag- 
netization allows to extend the same type of analysis to 
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Figure 2: The charge parity order parameter measured at half 
chain 0^p\L/2) as a function of the local interaction U . We 
have considered PBC's and finite chain lengths from L = 10 
to L = 50 in step of 4, a sequence for which the ground state 
of the Hubbard model is unique. We have plotted the curves 
for L = 10 (up triangles), L = 14 (circles), L = 18 (down 
triangles), L = 26 (rhombs). The filled squares represent the 
finite size scaling values to L = oo obtained by the fits shown 
in the inset. 



the LE phase, which is then characterized at any filling 
by bound pairs of single electrons with opposite spin. 

Based on the above scenario, we construct an approxi- 
mation scheme that aims at isolating the relevant degrees 
of freedom (charges) to describe the actual role of Op in 
the Hubbard model. Since the operator e*'^"^ — (— 1)"^ 
changes sign whenever the site j is singly occupied, no 
matter its spin orientation, we choose to represent the 
original electronic creation operators in terms of a 
spinless fermion f\ and Pauli operators uf, a = x,y,z 
acting on a spin part. The mapping, schematized in Ta- 
ble [T] is identified by the unitary transformation 



with Pj* — . Interestingly, the interaction term for 
the c-fermions simply becomes a chemical potential shift 
for /-fermions, namely ?7X]i'^iT"4 = ~ Li'^f)/^, 

where N — ^ rii^. According to this picture, the spin 

and pseudospin operators are sj'*' = fjfjCj and s'j^^ = 
fjfjtyj', conversely, we have (Tj = + J^. 



spinful fermion 
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It) 


U) 


Iti) 


spinless fermion ® tr-spin 


10)1+) 


ll)l+> 


|i>l-> 


|o)l-) 



Table I: Mapping from electrons to spinless fermions and a- 
spins. 



After the mapping the model in Eq.([2| becomes 
^ = - E [///^^'^ - 2(-l)7,/ /ji?., + H.c. 




(5) 



where Qij = {(Ji • cr^ + l)/2 is just the swap operator in 
the cr-spin state and P^^ = {l — (Ti-(Tj)/A is the projector 
onto the singlet. Notice that ([s]) is invariant under global 
CT-spin rotations. 

The form ([5| for the Hubbard model holds in arbitrary 
dimension, and its terms are quadratic with respect to /- 
fermions. Since Op ^ can be entirely expressed in terms 
of /i, a possible strategy consists on tracing out the a- 
spins by some mean-field approximation. In fact, ex- 
ploiting the symmetries of the Hubbard model one can 
easily realize that {Qij) = 1/2 is an exact identity on the 
states on which the hopping term in (|5| is non- vanishing. 



Moreover, we set the parameter a 



{Rij) in a phe- 



nomenological way by equating the ground state (GS) 
energy obtained from the spinless quadratic model with 
the exact energy coming from the Bethe-Ansatz solution 
jll| . Within this approximation Eq.(|5| is diagonalized in 
Fourier space, obtaining 



U{2N - L) 



cos k + ^JlQa'^ cos^ k + [7^/4 and 



keBZ 

with spectrum A^ 
rjk are the new fermionic modes. In the thermodynamical 
limit (TL), the energy density ecs at half- filling = 1 is 
given by ecs = ¥ - 2¥ ^-^2 Vl6a2 cos^ fc + C/74. It 
is interesting to observe that the model is gapless only for 
U = {), where for a = 1 6^5 assumes the exact value of 
the non-interacting case. For U > Q the number of singly 
occupied states Vf i& increasing and the pair-singlet states 
start to interact. 

We are interested in calculating the parity operator 
0^p\r) = (e'''5:}l-(«|-i))^ that can be rewritten as 



0^^\r) = lX{{2f]f,-l)) = lY{A,B, 



I i-\-r 



having defined = {f]+fj) and = ifj-f])- Making 
use of the Wick theorem. Op' (r) can be expressed as a 
determinant [12 1 
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where G is a block Toeplitz matrix of dimension (r + 1) x 
(r + 1), whose entries are the one-body correlation func- 
tions Gr = ((/j - fj){f]+r + fj+r))i whosc cxprcssions 
in the TL are 



U rTT cos(fcr) 

JO ^16a2 cos2(fe) + ;72/4' 

/ -|^\j4Q r^r cos(fcr) cos(fc) 

/ ' Jo ^16a2cos2(fc) + C/2/4' 



r even 



r odd 



with the property that Gq = 2i/f — l, Gr — 2Re((/J/j+r)) 

for r even and Gr ~ 2Re((/J/J_,_^)) for r odd. The blocks 
in ^ are of size 2x2. We must distinguish the cases 
of r even or odd, since they give rise to two different 
sequences. In particular, here we stick to the case r odd, 
where the block matrix is of even dimension. 

(c) 

The analytical calculation of Op in the TL for some 
values of U in the /-fermion approximation yields to the 
curve plotted in FigjS] evidentiating the expected non 
vanishing of the charge parity order for U > 0. The pa- 
rameter a has been determined by requiring ecsiU, a) ~ 
Sex{U), where Cex is the exact result Remarkably, 
such equality admits a solution for every U, which be- 
longs to a narrow interval below a = 1, as shown in 
the inset of FigjS] This means that in the pair-creation 
processes in Q the cr-spin state is very close to the sin- 
glet. In the limit U ^ 1 the energy becomes ecs{<^) ~ 
— 4a^/C/ that gives a{U oo) = \/log 2 w 0.83, by com- 
parison with the energy density of the Heisenberg model 
coming from the large- C/ expansion of the Hubbard model 

(c) 

a.t 1/ ~ 1. The result for Op is also quantitatively in ac- 
cordance with the DMRG data in the large U region, 
where our assumptions on the AF nature of short-ranged 
correlations |T3] are more justified. 

In conclusion, our work unveils that two (charge and 
spin) hidden parity string correlators play the role of or- 
der parameters for the gapped phases of the Hubbard 
model. In the bosonization approach these are found to 
be asymptotically finite only in the corresponding gapped 
MI and LE phases, and vanish with the gap at BKT tran- 
sition point. The result is cleanly confirmed by DMRG 
numerical analysis. The emerging physical scenario is 
that of an insulator in which bound pairs of doublons and 
holons move in a AF background of single electrons. In 
the LE regime the role of doublons and holons and that 
of up and down electrons are exchanged. The picture 
allows to derive an effective free spinless fermion model 
which captures correctly the presence of non local order, 
and its vanishing at the transition. The spinless model 
is exact in the limit of large U thus complementing the 
standard strong-coupling description with t — J model. 

The parity order is suitable for experimental detection 
by high resolution imaging [B] in ultracold Fermi gases. 
Possibly, some of the features described here could per- 
sist in two dimensions, where the localization of bound 
pairs could take place along one dimensional stripes. The 
scenario is quite suggestive also from the perspective of 



0.6 

o 

0.4 









d' / 0.95 
□ ' J 

- / f a„ 









u 



10 12 14 



Figure 3: The parity order Op' calculated according to Eq.(|6| 
as a function of U (dashed line) in the MI phase of the effective 
model Q with a as in the inset figure, determined by tuning 
the spinless fermion energy to the exact Hubbard value. The 
results are compared with the numerical curve (continuous 
line) obtained for the Hubbard model and shown in Fig|2] 



high- Tc materials: the presence of bound doublon-hole 
pairs in the undoped insulator could play a role upon 
doping in the onset of the superconducting phase. 

The present analysis could be further exploited to ex- 
tended Hubbard models |14| IT5| , to describe other topo- 
logically ordered phases; noticeably, the fully gapped 
phase characterized by non vanishing charge and spin 
gaps should correspond to the non vanishing of both 
oP's. Work is in progress along these lines. 

M.R. acknowledges support from the EU-ERC project 
no. 267915 (OPTINF) . 
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